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ABSTRACT 

In this work we simulate the 50 — 200 MHz radio sky that is constrained 
in the field of view (5° radius) of the 21 Centimeter Array (21CMA), a low- 
frequency radio interferometric array constructed in the remote area of Xinjiang, 
China, by carrying out Monte-Carlo simulations to model the strong contami- 
nating foreground of the redshifted cosmological reionization signals, including 
emissions from our Galaxy, galaxy clusters, and extragalactic discrete sources 
(i.e., star- forming galaxies, radio-quiet AGNs, and radio- loud AGNs). As an 
improvement of previous works, we consider in detail not only random varia- 
tions of morphological and spectroscopic parameters within the ranges allowed 
by multi-band observations, but also evolution of radio halos in galaxy clusters, 
assuming that relativistic electrons are re-accelerated in the intra-cluster medium 
(ICM) in merger events and lose energy via both synchrotron emission and in- 
verse Compton scattering with cosmic microwave background (CMB) photons. 
By introducing a new approach designed on the basis of independent component 
analysis (ICA) and wavelet detection algorithm, we prove that, with a cumu- 
lative observation of one month with the 21CMA array, about 80% of galaxy 
clusters (37 out of 48 clusters assuming a mean magnetic field of B = 2 fiG in 
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the ICM, or 15 out of 18 clusters assuming B = 0.2 /xG) with central brightness 
temperatures of > 10 K at 65 MHz can be safely identified and separated from 
the overwhelmingly bright foreground. By examining the brightness temperature 
images and spectra extracted from these identified clusters, we find that the mor- 
phological and spectroscopic distortions are extremely small as compared to the 
input simulated clusters, and the reduced x 2 °f brightness temperature profiles 
and spectra are controlled to be < 0.5 and < 1.3, respectively. These results 
robustly indicate that in the near future a sample of dozens of bright galaxy 
clusters will be disentangled from the foreground in 21CMA observations, the 
study of which will greatly improve our knowledge about cluster merger rates, 
electron acceleration mechanisms in cluster radio halos, and magnetic field in the 
ICM. 

Subject headings: early universe — galaxies: clusters: general — radio contin- 
uum: general — reionization — techniques: image processing 

1. INTRODUCTION 

Galaxy clusters are the largest virialized systems in the Universe. Each of them contains 
about 10 14 — 10 15 M matter in the form of galaxies, hot gas and dark matter, and all 
these components have been studied in depth in optical, X-ray, and high-frequency (> 1 
GHz) radio bands (e.g., Sarazin & Kempner 2000; Rosati et al. 2002; Popesso et al. 2004; 
Ferrari et al. 2008). However, the physical conditions of magnetic fields and high energy 
particles, which permeate the extremely hot intra-cluster medium (ICM) and have been 
playing important roles in the evolution of the ICM (e.g., Schekochihin et al. 2005; Lazarian 
2006), are still poorly understood. In the study of giant radio halos, which are always found 
in central regions of galaxy clusters and usually span on Mpc scales, these lacked knowledges 
actually become crucial, because most of the emission of the radio halos can be ascribed to 
synchrotron radiation from relativistic electrons moving in the /iG magnetic fields during or 
after a merger event (e.g., Ferrari et al. 2008; Brunetti et al. 2008). A variety of electron 
acceleration mechanisms and magnetic field distributions (e.g., Enfilin & Gopal-Krishna 
2001; Brunetti et al. 2001, 2004; Cassano & Brunetti 2005) have been proposed to construct 
dynamic models to describe the behavior of relativistic electrons, and these models are to 
be testified by the upcoming facilities working in the low-frequency (< 1 GHz) radio band, 
such as the 21 Centimeter Array (21CMA; Peterson et al. 2005), the Low Frequency Array 
(LOFAR; Rottgering et al. 2003), and the Murchison Wide-field Array (MWA; Morales & 
Hewitt 2004). 
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As of today, only < 50 detections of radio halos in the GHz radio band have been re- 
ported (e.g., Cassano et al. 2006; Giovannini et al. 2009; Brunetti et al. 2009). However, 
since the characteristic decay time of synchrotron emission becomes significantly longer when 
the frequency goes lower (for example, the synchrotron decay time at 150 MHz is approx- 
imately three times longer than that at 1.4 GHz), much more radio halos are expected at 
frequencies of a few hundred MHz. For example, Cassano et al. (2010) estimated that more 
than 350 giant radio halos will be detected at 120 MHz by LOFAR in northern hemisphere 
and at Galactic latitudes |6| > 10°. And this expectation has been partially confirmed by 
recent Giant Metrewave Radio Telescope (GMRT) observations at 610 MHz (e.g., Venturi 
et al. 2007, 2008). Clearly, tight observational constrains on the formation and evolution of 
radio halos of galaxy clusters will be obtained in the next decade. 

In most previous works on future low-frequency radio researches, galaxy clusters were 
only treated as part of the useless contaminating foreground in front of the redshifted 21 
cm emission signals that come from the cosmic reionization era, and were simply filtered 
out, together with the foreground emissions from our Galaxy and extragalactic discrete 
sources (e.g., Wang et al. 2006; Jelic et al. 2008). Little effort has been made to study the 
feasibility of separating the cluster component from the luminous foreground. Meanwhile, 
we find that in the existing numerical works models constructed for galaxy clusters are often 
oversimplified by either introducing a universal spectral index for the radio flux density of 
all cluster halos (a CL = 1.2; e.g., Gleser et al. 2008; Jelic et al. 2008), or, sometimes, by 
assuming a flat profile (e.g., Jelic et al. 2008) for the spatial distributions of the radio surface 
brightness. Similar oversimplifications are often found in the simulations of the emissions 
from our Galaxy and extragalactic discrete sources, too. 

In order to construct the contaminating foreground in front of the redshifted cosmic 
reionization signals with reasonably higher spatial and spectral accuracies, by fully taking 
into account the random variations of model parameters in the ranges allowed by the obser- 
vations, and to investigate whether or not the cluster component (or, at least, brightest clus- 
ters) can be successfully disentangled from the bright foreground, we carry out Monte-Carlo 
simulations to model the Galactic, cluster, discrete source, and the reionization components, 
with which we create the 50 — 200 MHz radio sky with high fidelity as possible as we can. 
Then we attempt to separate galaxy clusters and extract brightness temperature profiles and 
spectra from the brightest clusters, whose central brightness temperature are brighter than 
10 K at 65 MHz, by introducing an approach designed on the basis of the independent com- 
ponent analysis (ICA; Comon 1994) and wavelet detection algorithm (Shensa 1992). All the 
simulations are performed for the 21CMA array that covers a field of view (FOV) of about 
5° — 10° around the North Celestial Pole (NCP), depending on the frequency. Throughout 
the work we adopt H =71h n km s" 1 Mpc^ 1 , fi M = 0.27, fi A = 0.73, and tt b = 0.044 (e.g., 
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Spergel et al. 2003). 

2. SIMULATION OF THE LOW-FREQUENCY SKY 

In this section, we simulate the contributions that come from our Galaxy, galaxy clusters, 
extragalactic discrete sources, cosmic reionization regions, and antenna array, respectively. 
All the simulated results are plotted in a circular FOV with a 5° radius around the NCP, 
as will be observed with the 21CMA array, although in our simulations we firstly create 
corresponding images with a radius of 6°, so that we are able to correct the boundary 
effect. We project the sky field to a flat plane using the azimuthal equidistant projection, 
and the obtained images are overlaid with grids of 1024 x 1024 pixels, each pixel covering 
approximately a 0.6' x 0.6' patch. 

2.1. Our Galaxy 

In the low-frequency radio band the emission of our Galaxy is luminous and consists 
of synchrotron, free-free, and dust components (e.g., Jones et al. 1998a; Platania et al. 
2003; Gleser et al. 2008). Of these, synchrotron emission is the dominating component, 
and is responsible for ~ 98.5% of the total Galactic contamination in 100 — 200 MHz. Free- 
free emission from the diffuse ionized hydrogen in the interstellar medium (ISM) is the next 
important contaminating source that contributes up to ~ 1.5% of the total Galactic emission 
in 100 — 200 MHz. Dust emission is mainly restricted on the Galactic plane and is usually 
negligible in 100 - 200 MHz (Shaver et al. 1999; Platania et al. 2003), thus we don't take it 
into account in this work. 

The diffuse Galactic synchrotron emission exhibits visible substructures in the total 
(Stokes I) intensity map, which can be attributed to the inhomogeneous spatial distributions 
of the magnetic field (e.g., Gaensler et al. 2001) and the ISM (e.g., Wieringa et al. 1993; Gray 
et al. 1999; Haverkorn et al. 2000). As of today, the smallest substructures were resolved 
on the ~ V scales in the 1.4 GHz Southern Galactic Plane Survey (SGPS; Tucci et al. 
2002), whereas in the low-resolution Westerbork Telescope (WSRT) observations lO'-scaled 
substructures at 150 MHz were reported as a common feature (Bernardi et al. 2009). In 
our simulation, we construct the two-dimensional Galactic synchrotron map from observation 
results (Giardino et al. 2002) directly, rather than the modeling of three-dimensional Galactic 
structures (e.g., Jelic et al. 2008; Sun et al. 2008; Sun & Reich 2009). By adopting the 
21CMA FOV corresponding region of the brightness temperature image obtained at 408 
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MHz (T 4 G s 8 y 

MHz( r )' where r is the two-dimensional position; Giardino et al. 2002), on which 
substructures with scales larger than 0.85° are reproduced based on the 408 MHz all sky 
map (Haslam et al. 1982) and substructures with smaller scales (but > 2') are simulated 
using the extrapolation of the Galactic emission power-spectrum, we obtain the simulated 
Galactic synchrotron emission images 

T& yn (r) ( w here v is the frequency; see an example 

T 6?MHz( r ) in Fi g- la ) as 



where a Gsyn (r) ~ 2.7 — 2.9 is the corresponding spatially varying synchrotron temperature 
spectral indeXLj (Fig- 16) calculated from the data of the 408 MHz and 1.4 GHz surveys 
(Giardino et al. 2002). Thus we obtain the two-dimensional Galactic synchrotron emission 
distribution T„ syn (r), which ranges from mean value of 5213 K (root mean square, RMS 
hereafter, is 712 K) at 50 MHz to 106 K (RMS is 10 K) at 200 MHz. 

The Galactic free-free emission component cannot be directly observed with today's 
technology, because the Galactic synchrotron component dominates the emission at all fre- 
quencies lower than 10 GHz, meanwhile the dust thermal emission becomes overwhelming at 
higher frequencies. However, hydrogen recombination lines (Ha line emission in particular) 
provide a novel way to trace the free-free emission, since they originate in nearly the same 
region (Smoot 1998; Reynolds & Haffner 2000). In fact, the brightness temperature distri- 
bution of the Galactic free- free emission component at 30 GHz TgjJ GHz (r) has been related 
to the Ha intensity distribution Iff a (r) (Finkbeiner 2003) in the form of 

where 1 Rayliegh = 10 6 /47r photons s _1 cm -2 sr -1 (Reynolds & Haffner 2000). Using this 
relation, we derive the T^ s GRz (r) map of the 21CMA FOV from the corresponding region 
of Iho{^) survey image (Finkbeiner 2003). By assuming that the brightness temperature 
spectrum of the Galactic free- free emission follows a broken power-law shape, i.e., T^ s oc 
u~ a T B , where a% s = 2.10 at v < 10 GHz (Shaver 1999) or a% s = 2.15 at v > 10 GHz 
(Bennett et al. 2003) is the free-free temperature spectral index, we theoretically obtain the 
two-dimensional Galactic free- free emission distribution T^ G (r) (see an example T 6 ^ ff MHz (r) 
in Fig. lc), which ranges from mean value of 9 K (RMS is 2 K) at 50 MHz to 0.5 K (RMS 
is 0.1 K) at 200 MHz. 



1 According to Ray leigh- Jeans approximation, temperature spectral index ax is related to the flux density 
spectral index a via ar — a + 2. 
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2.2. Galaxy Clusters 

2.2.1. Mass Function and Three-dimensional Spatial Distributions 

The mass distribution of dark matter halos at a given redshift z can be derived in 
the frame of the Press-Schechter model (P-S model; Press & Schechter 1974), which has 
been extensively used in a variety of works (e.g., Zhang & Wu 2003; Zentner et al. 2005; 
Fedeli et al. 2006) since it was proposed. According to the P-S model, the initial per- 
turbation of dark matter density at a given position R, i.e., p(R), is characterized by the 
density contrast 5(R) = ^ with its probability p(S) following the Gaussian distribution 

= y/2mr{M) e~ 52 1 2a2 where a(M) is the variance in the density field and M is the 
mass of dark halo. For a smoothed density field with a constant smoothing scale of R sm , 
inside which the total mass is M, the variance in the density field is o~(M) = <((5|, sm (R)) = 
2^ /^fl sm (^) k 2 dk, where k is the wave number, WR sm {k) = e~ k R ^l 2 is the fil- 
ter function for Gaussian distribution, and P(k) = \5k\ 2 is the primordial power-spectrum. 
Based on these the dark halo mass function can be expressed as 

where S c (z) = is the halo collapse time, D(z) is the normalized growth factor, and 
S c (0) = 0.15 (12tt) 2 / 3 fi^ 0055 , in which fi M is the total matter density. 

Assuming a mass fraction of 83% for dark matter (e.g., Drexler 2007) and a lower mass 
limit of 2 x 10 14 M (e.g., Holder et al. 2001; Voit 2005) for galaxy clusters, we calculate 
the number density distribution of galaxy clusters as a function of mass and redshift as 
n CL _ n (j\^fCi^ z ^ (M CL is the cluster mass) according to the P-S model, and the redshift 
distribution of galaxy clusters as n CL = n CL (z) (Fig. 2). By examining the derived n CL (z), 
we find that galaxy clusters mainly distribute in the redshift range of < z < 3 (e.g., Holder 
et al. 2001), which agrees with the observations, and 1084 cluster-level halos are expected in 
the 21CMA FOV. In the work that follows we assign random M CL and z to these 1084 clusters 
according to the n(M CL , z) distribution derived above, and determine the coordinates of the 
clusters on the sky plane (i.e., R.A. and Dec) according to the predictions of the theoretical 
cluster power spectrum in redshift range of < z < 3 (Majumdar & Mohr 2004). 



2.2.2. Spectral Model 



Galaxy clusters that contain giant, diffuse radio halos are observed to exhibit a strong 
correlation between their X-ray and radio powers (e.g., EnBlin & Rottgering 2002), whereas 
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those contain no or weak radio halos are found to disobey this relation, showing a significantly 
lower L radio /L x ratio (e.g., Cassano 2009), where L r ^ Q and L x are the radio and X-ray lu- 
minosities for the clusters, respectively. According to the standard turbulent re-acceleration 
model, the primordial relativistic electron population accelerated in a merger event follows 
a power-law energy distribution N(pj, t)\t=o oc 7(5~ r , where the Lorentz factor is 7 ^> 1000, 
70 = 7|t=o, an d the electron energy spectral index is T > 3 (e.g., Ferrari et al. 2008). These 
electrons lose energy not only via synchrotron emission when moving in the magnetic field 
(B), but also via inverse Compton scattering through interacting with the cosmic microwave 
background (CMB) photons. The latter process can be described in a similar way to the 
synchrotron emission if an equivalent magnetic field of -Bcmb = 3.2(1 + z) 2 /xG is introduced 
(e.g., Parma et al. 1999), thus the electrons lose their energy on a time scale of 



(Ferrari et al. 2008). Taking z = 0.2 (around which most of the giant radio halos are found) 
and the typical magnetic field strength B = 2 /xG (e.g., Kim et al. 1990) or 0.2 /xG (e.g., 
Ferrari et al. 2008; Rephaeli & Gruber 2003; Rephaeli et al. 2006) for example, we find that 
the emission of the radio halo at 1.4 GHz will fade away after r ~ 0.06 Gyr or ~ 0.02 Gyr for 
B = 2 /xG or 0.2 /xG, respectively. At 100 MHz, the corresponding decay time will increase 
significantly to r ~ 0.2 Gyr or ~ 0.1 Gyr for B = 2 /xG or 0.2 /xG, respectively. This means 
that the spectrum of a radio halo will evolve into a soft state quickly, and finally fade away 
in the interval between two merger events (typically 2.7 Gyr; Mathiesen & Evrard 2001). 
As a result, typical radio temperature spectral indices of galaxy clusters are a^ L u > 3 (Fig. 
3 a), which are significantly larger than those of other foreground components, including our 
Galaxy (§2.1) and extragalactic discrete sources (§2.3). In this work, we assign a random 
evolving time t E between and 2.7 Gyr for each cluster, and calculate the low- frequency 
flux density F^ L (t E ) from its power density P^ L (t E ) (see details in Appendix) by assuming 
that each cluster owns a radio halo produced in the last merger and the initial 1.4 GHz 
power density -P^ghz^ = 0) obey the observed correlation between the 1.4 GHz power 
density and X-ray temperature Tx- In order to derive the approximation of PfJ GHz — Tx 
relation, using the data from Cassano et al. (2006), we repeatedly fit the observed data with 
maximum likelihood method and exclude the clusters fall out of the 3a lower limit, as is 
shown in Figure 3b. Then we have 
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ergs s 1 Hz 
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(5) 



As a comparison, on the figure we also superpose the data of three new cluster samples stud- 
ied by Giovannini et al. (2009), and find that all of them locate below the line described by 
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equation (5). The X-ray temperatures of the clusters are deduced from the mass-temperature 
relation provided by Arnaud et al. (2005), i.e., 

CL 1 1/1-72 



T x = 5 



h(z)M' 



keV, 



(6) 



5.34 x 10 14 M Q 

where h(z) = a/^m(1 + z) 3 + 1 — Qm is the Hubble constant normalized to its local value. 



2.2.3. Morphologies 



Since for galaxy clusters there exists a linear correlation between the radio and X-ray 
surface brightness distributions (e.g., Govoni et al. 2001), we simulate the radio surface 
brightness density distribution by applying the (3 model (Cavaliere & Fusco-Femiano 1976) 



CL 



1 



-3/3+0.5 



(7) 



where r c is the core radius, /3 is the index, and Sq L (u) is the central surface brightness 
density. The core radius r c is derived from the cluster's X-ray luminosity in 0.5 — 2.0 keV 

(£x,o.5-2.okcv) by 

176 / £x,0.5-2.0kcV 



r r . = 



0.2 



kpc 



(8) 



h 71 V.5 x 10 44 ergs s" 1 

(Jones et al. 1998b), and £x,o.5-2.okeV is calculated from the bolometric (0 — 20 keV) lumi 
nosity using 

T, 



L Xibolo = 12.44 ± 1.08 x 10 



44 



'6 keV 



x ^2.64±o.27 ^ _|_ z y->: -\ 



,+0.26 

-o.27 er g S s 



(9) 



(Lumb et al. 2004), assuming the metal abundance as a function of redshift as given in 
Ettori (2005). On the other hand, we calculate the parameter f3 by applying equation (2) of 
Evrard et al. (1996) 

- 1 M CL 



P = 8.85 x 10 



-is 1 + (Wr c ) 2 ( T x 



(r 2 oo/r c ) 



keV 



r 200 

Mpc 



M ' 



(10) 



where r 2 oo is the virial radius within which the mean interior over-density is 200 times of the 
cosmic critical density at the redshift considered. When r c and /3 are determined, Sq L u can 
be constrained from cluster's flux density F^ L (ts), which is given in §2.2.2, by 



/T200 
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/ ^rS$% 
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-3/3+0.5 



dr. 
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Thus, having derived the parameters r 20 o, coordinates, r c , (3, and Sq L u of each cluster, 
we generate the two-dimensional surface brightness density map of each galaxy cluster at v 
using equation (7), overlay all the maps together to construct the simulated map of cluster 
component in the 21CMA FOV, and convert the map into simulated cluster brightness 
temperature image Tg^(r) (see examples at 65 MHz in Fig. 4) with the Rayleigh- Jeans 
approximation 

T B>)=(^ B )'V(r). (12) 

2.3. Extragalactic Discrete Sources 

In most simulations on the foreground of the redshifted 21 cm reionization signals (e.g., 
Wang et al. 2006; Santos et al. 2005; Jelic et al. 2008), it is assumed that all the extragalactic 
discrete sources are point-like and possess the same radio spectrum in the form of F u oc z/~ ar , 
whereas the observations indicate a more complicated situation. Following the works of 
Wilman et al. (2008; W08 hereafter) and Snellen et al. (2000; S00 hereafter), we classify 
the extragalactic sources into four types, i.e., 1) star-forming galaxies, including normal 
star-forming galaxies and starburst galaxies; 2) radio-quiet AGNs; 3) Fanaroff- Riley class I 
(FRI) and class II (FRII) AGNs; and 4) GHz peaked spectrum (GPS) and compact steep 
spectrum (CSS) AGNs. We simulate the flux densities, spatial structures, spectra, and 
angular clusterings of types (l)-(3) sources following Wilman et al. (2008) and references 
therein, and those of type (4) sources following Snellen et al. (2000) and reference therein. 
To be specific, we start with the flux density functions that were obtained for types (l)-(3) 
sources at a fiducial frequency Uf = 1.4 GHz (Fig. 4 of W08), and the evolving luminosity 
function for type (4) sources at Vf = 5 GHz (Eq. (21) of S00), respectively, to derive the flux 
densities of the sources that should appear in the 21 CM A FOV, and simulate their spatial 
structures to create the brightness temperature images at Vf. Then with the information of 
their spectral properties, we simulate the brightness temperature images in 50 — 200 MHz 
band for each type of extragalactic discrete sources (see examples at 65 MHz in Fig. 5), and 
obtain the brightness temperature images of the extragalactic discrete component Tg ls ^(r). 
The simulation details in our process are as follows. 

Star-forming galaxies 

(i) Each source is assumed to be disk-like and located on the sky plane with a random 
orientation. For a normal star-forming galaxy, the linear size D is related to the 1.4 GHz 
luminosity by log(^) = 0.221og( w £ I ^_P^_ 1 ) - log(l + z) - 3.32(±0.18) (Eqs. (7)-(9) of 
W08). For a starburst galaxy, we assume that D — (1 + z) 2 - 5 kpc (Eq. (10) of W08) out to 
z = 1.5 and D = 10 kpc for z > 1.5. 



- 10 - 



(ii) The emission of each star-forming galaxy consists of one thermal free-free and one 
non-thermal synchrotron component, which are constrained by the ratio of the total lumi- 
nosity to the luminosity of the thermal component of 1 + 10 (q^) °' 65±0 ' 10 (W08). 

(iii) The spectrum of each source is modeled with two power-laws as L v oc (q^) ° 10 + 
10 (gSz) ° 75 ° 10 (^) °^ W08), where the two terms represent the thermal free-free and 
non-thermal synchrotron components, respectively. 

Radio-quiet AGNs 

(i) Sources are all observationally point-like. 

(ii) The spectrum of each source can be modeled with a power-law model as F v oc z/~ ar 
(W08) that is normalized by using the 1.4 GHz luminosity. 

FRI & FRII AGNs 

(i) Each source is assigned with a linear size that follows a random, uniform distribution 
[0, Dq(1 + z) -1 ' 4 ] at a given redshift, where D = 1 Mpc. Each source possesses a point-like 
compact core, and two extended lobes that can be modeled as ellipses with the axial ratio 
drawn from a uniform distribution [0.2, 1] (W08). 

(ii) According to the orientation-based unification model, the ratio of the 1.4 GHz flux 
density of the core to that of the lobes is F^ e GH jF^ e GHz = R v{ B(6), where .Rrf = 10 x is 
the rest-frame ratio with x following a Gaussian distribution N(x, 0.5), and B{6) = 1/2[(1 — 
■sj (7 2 — l)/7cos#)~ 2 + (1 + a/(7 2 — l)/7cos#)~ 2 ] with 7 being the Lorentz factor of the 
jet and 9 being the angle between jet axis and line-of-sight (i.e., inclination angle) that 
distributes uniformly between and n. We set (x, 7) = (—2.6,6) and (—2.8,8) for FRI and 
FRII sources, respectively (W08). 

(iii) The spectrum of the core emission is estimated as logF^ ore = ao + ailog(^^) + 
a 2 log 2 (^^), with ai = 0.07 ± 0.08, a 2 = —0.29 ± 0.06 and a being calibrated by applying 
F^Xghz- The spectrum of the extended lobe emission is parameterized as i^ obc oc z/-°- 75 
(W08). 

GPS & CSS AGNs 

(i) Sources are all observationally point-like. 

(ii) The spectral shape is found to have a turnover below several GHz, which can be 
described as 



it.GPS,CSS 
^GPS,CSS _ v v ( V 



J± f-y* i- e -(%>- « (13) 

1 — e~ L Up L J 

(Eq. (1) of Snellen et al. 1998), where t k = 0.51±0.03 and t n = -0.73±0.06 are the optically 
thick and thin spectrum indices (de Vries et al. 1997), respectively, and v v is the observed 
turnover frequency. The turnover frequency varies with log Vp ^ z) = -0.21(±0.05) - 
0.65(±0.05)log(^) (Eq. (4) of W08), in which the projected linear size of the source D is 
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independent from redshift but follows the scale function between D and the source number 
N, 

JAr/Jn ( D- 1 0.01 kpc < D < 0.355/^ kpc M 

(Fanti et al. 2001; O'Dea 1998; Stawarz et al. 2008). 

The total flux density functions of our simulated extragalactic discrete sources are shown 
in Figure 6. 



2.4. Effects of Instrument 

The 21CMA array is a low-frequency radio interferometer array constructed in a remote 
area of Xinjiang, China. The array consists of 81 antenna pods that are distributed in a 
"T" shape with the north-south and east-west baselines are 6 km and 4 km, respectively, 
and each pod is an assembly of 127 logarithmic periodic antennas. The array is operated in 
50-200 MHz with the minimum bandwidth of 49 kHz. 

In radio interferometric observations, the systemic uncertainty for one image pixel can 
be expressed in terms of brightness temperature as 

A 2 

AcfffibeamV Avm(n ~ 1) 

(Thompson et al. 2001), where T sys is the system temperature, A is the wave length, A e s is the 
effective area of one antenna pod, f2b eam is the beam solid angle of the interferometric array, 
Av is the bandwidth of one frequency channel, r is the integral time, and n is the number of 
antenna pods. The 21CMA beam varies with the wavelength A, e.g., the full width at half 
maximum (FWHM) of the beam is about 3.2' (corresponding to fibcam = 6.8 x 10~ 7 sr) at 65 
MHz, and about 1.6' (n bcam = 1.7 x 10 -7 sr) at 150 MHz. For 21CMA, we have T sys = 300 K, 
A cS = 216 m 2 , n = 81, thus, 

_ f A A 2 (\ MHz\ 1/2 /30 days\ 1/2 fl.7 x lO" 7 srA 

Supposing a typical survey of r = 30 days and Av = 1 MHz, for each frequency channel we 
generate the instrument noise as white noise and add it into the simulated 21CMA map. 

By mixing the simulated brightness temperature images of all celestial and noise com- 
ponents, we create the final sky maps constrained in the 21CMA FOV (Tj^^r), where 50 
MHz < v < 200 MHz, Av = 1 MHz), assuming that in the ICM of galaxy clusters there 
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exists a uniform magnetic field B = 2 /iG and 0.2 /iG, respectively. We apply a universal 
Gaussian kernel that is required by following component separation approach (see details 
in §3), whose FWHM= 3.2' corresponds to the beam solid angle fibeam = 6.8 x 1CT 7 sr at 
65 MHz, to approximate the point spread function (PSF) of the observation, and convolve 
the simulated sky maps that have been cumulated within selected frequency channels (cen- 
tral frequency v = 61, 62, 63, 189 MHz, respectively) to have a channel bandwidth of 
Av = 1 MHz for each. These fake 21CMA survey maps are marked as 7^ ot 1 f 1 (r) hereafter 
(see examples with the central frequency of v = 65 MHz in Fig. 7a k, 7b). 

For the convolved sky maps, we calculate the confusion level with three different meth- 
ods, and have crosschecked the results to confirm if they are consistent with each other (Table 
2). In the first method, the confusion level is defined in a traditional way in terms of source 
flux density, over which there is 0.1 source per synthesized beam on average (e.g., Rottgering 
et al. 2006). We adopt the source flux density function (N(> F u )) that derived from our 
simulation (§2.3 and Fig. 6) and the beam (FWHM=3.2') of our simulated survey maps to 
calculate the confusion level at 65 MHz and 150 MHz. In an alternative way, we also replace 
our source flux density function with the literature source flux density function that has been 
used to evaluate the confusion limit of LOFAR (Appendix B and Figure 4 of Rottgering et 
al. 2008). In the second method, we employ the photometric criterion for confusion level 
as mentioned in Dole et al. (2003). From their equations (5)- (6) and related description, 
we deduce the confusion level to be [(3 + a)/(r2 b cam-^o? 2 )] 1 ^ a+1 \ where a and F are the 
index and normalization of the differential flux density function (—dN/dF u ), respectively, 
and q = 3 is the minimum signal to noise ratio for source detection. In the third method, we 
directly detect the sources on our simulated survey maps with the tool SExtractor (Bertin & 
Arnouts 1996) at a 3a confidence level. We derive the confusion level by counting the source 
number as a function of flux density, which is the approach employed by Vaisanen et al. 
(2001). Meanwhile, we also calculate the confusion level with the instrument beam size of 
the 21CMA array at 65 MHz (FWHM=3.2' also) and 150 MHz (FWHM=1.6'), respectively, 
by repeating the steps above. Since all the calculations above are at a 3cr confidence level, 
the value of confusion noise should be 1/4 of the confusion level. 



3. IDENTIFICATION OF BRIGHT GALAXY CLUSTERS 

In this section, after extremely bright discrete sources are masked, we apply the ICA 
separation and wavelet detection algorithms to these maps to identify and study bright 
galaxy clusters, as described below. 
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3.1. Step 1: Masking Extremely Bright Discrete Sources 

According to our simulations, except for our Galaxy, extragalactic discrete sources con- 
tribute most of the foreground in 50 — 200 MHz in the 21CMA FOV. Of these, about 
60 sources are extremely bright, and should be discarded before identifying and analyzing 
weaker components. Here we introduce a masking method, which has been widely used in 
the CMB separation works (e.g., Bennett et al. 2003; Leach et al. 2008), to exclude map 
pixels that contain extremely bright discrete sources. To be specific, we adopt a brightness 
temperature criterion of T^ ot ^ l (r) > 1.5 x 10 4 K in all selected frequency channels to exclude 
circular regions with a radius of one FWHM (3.2') centered at all extremely bright sources. 
In the analysis that follows, if a pixel is masked in any one frequency channel, it is also 
masked in all other channels. As a result, about 0.59% of the 21CMA FOV is masked. 



3.2. Step 2: ICA Separation and Wavelet Detection 

3.2.1. Raw Separation with FastICA 

As a blind signal separation method, ICA only assumes that the maps of original celestial 
and noise components are independent from each other in frequency space, and no more 
than one of them is Gaussian. These assumptions are applicable to our case. Thus, in terms 
of the ICA algorithm, the relation between the mixed maps obtained in different frequency 
channels, all of which have been whitened (Hyvarinen 1999) as a pre-procession, and original 
component maps can be expressed as 

x = As, (17) 

where x = {x±, x 2 , x m } T is the vector constructed by the mixed maps at to layers of 
frequency channels (i — 1, , 2 , to), s = {si, s 2 , s n } T is the vector constructed by n {n 
is the number of celestial components, n < to) layers of original celestial maps, and A is an 
to x n dimensional transition matrix. Since from equation (17) we have 

s = A _1 x, (18) 

where A" 1 is the inverse of A, the original celestial maps can be derived when A -1 is 
available. In order to determine A -1 , we introduce a random Gaussian field /x, with zero mean 
and unit variance, which has the same image size as the fake 21CMA maps, and maximize 
the difference between two fields /i and wfx, where Wj is the ith (i — 1, ,2 , n) row of 
the test matrix W with unit norm, by maximizing the approximation of nongaussianity, 

n 

Jam = J2i E i G ^M - E { G M]} 2 (19) 
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(Hyvarinen 1999), where G is the nonquadratic operator for log cosh, and the expectation 
E are calculated over the whole image. When the maximum is found, we have W = A -1 . 

As shown in equation (17), the ICA algorithm reads in and treats the input mixed maps 
as vectors, the elements of which are related to the pixels of the corresponding maps. The 
way to establish the mapping between a map and its corresponding vector (e.g., row- by-row, 
column-by-column, or even randomly) is arbitrary, but once determined for one map, it 
should keep the same for all the maps. Thus the spatial relations among pixels (i.e., image 
domain features) are broken and are not considered by the ICA algorithm. This is different 
from the traditional way to detect discrete sources on a map, for which the confusion noise 
is essentially important. 

ICA has been proved to be an effective separation technique in astronomical imaging 
analysis works, to reconstruct the sources confused by other components or even sunk into 
the confusion noise. In De Bernardis et al. (2009; the B-Pol proposal), the authors proposed 
to adopt ICA technique to extract dim point sources, which cannot be directly detected from 
the residual confusion noise but would dominate over the CMB B-mode signal. Baccigalupi 
et al. (2000) applied the ICA technique to a 15° x 15° sky patch at high Galactic latitude 
(I > 82.5°), and successfully reconstructed both the CMB signals with an RMS error of about 
1%, and the weaker (about two orders of magnitude fainter) Galactic components with an 
accuracy of ~ 10%. Maino et al. (2006) and Bottino et al. (2008) applied an optimized 
ICA (FastICA; Hyvarinen 1999) technique, which adopts fixed-point algorithm in ICA to 
maximizing Jq more efficiently, on the whole sky Wilkinson Microwave Anisotropy Probe 
(WMAP) map to extract the CMB signals and separate the multiple foreground components. 

In this work we attempt to use the FastICA package embedded in the Modular Toolkit 
for Data Processing (MDP; Zito et al. 2009) to disentangle the cluster component from 
the fake 21CMA survey maps. To reduce the complexity of the calculation, we divide the 
21CMA FOV into one circular and 24 pie patches (Fig. 7c & Id) that have the same 
areas. In order to correct the boundary effect, we enlarge each patch by at least 10' so 
that it slightly overlaps all the adjacent patches. We carry out separation calculation in 
each patch by feeding the fake 21CMA maps T^ otal (r) obtained in the v = 65, 80, 170, 
and 185 MHz channels into the FastICA program. In the calculation, we assume that 
each celestial source bears similar morphologies in different frequency channels, due to the 
source correlation coefficient between any two channels is > 0.9, which is required by the 
algorithm. We start the FastICA calculation in Patch I with an arbitrarily chosen test 
matrix Wj. When the maximum Jq is found as W™ ax (= Aj -1 ), we feed W™ ax into the 
FastICA calculation of Patch II as the original test matrix W n , because this will greatly 
reduce the iteration steps to find Wfj ax (= Afj 1 ). After repeating this progress for all patches 



- 15 - 



to find all W i max (= Af 1 ; % = I, II, III, XXV), we obtain raw separated images T^J(r) 
(Src=Gsyn, Gff, CL, and Disc for Galactic synchrotron, Galactic free-free, galaxy cluster, 
and discrete source components, respectively) using equation (18). 

In each frequency channel the output raw separated images for different components are 
rescaled to guarantee that the cumulative brightness temperature at each image pixel equals 
to the corresponding input value. To deal with masked regions containing extremely bright 
discrete sources, we fill each of them with the mean brightness temperature value calculated 
from the surrounding pixels after the raw separation is finished. 

3.2.2. Identification of Bright Galaxy Clusters 

In order to locate galaxy clusters on the raw separated images ^B L 65MHz'( r )> we introduce 
a detection algorithm based on the wavelet transform (WT; Slezak et al. 1990), which enable 
us to discriminate structures as a function of scale, especially those small-scale ones embedded 
within the larger ones. To be specific, we adopt the a Trous wavelet transform algorithm 
(Shensa 1992; Vikhlinin et al. 1997; Gu et al. 2009), which is a nonorthogonal version 
of discrete wavelet transform. We employ a Gaussian-shaped convolution mask h(cr, I) = 
e -i 2 /(2a 2 ) / /(2 7ro - 2 ) ) and define the corresponding wavelet kernel as W(l) = h(a,l) — h(2o-,l), 
where a is the scaling parameter and / = 2 n pixels (n = 1, 2, . . .) is the scale. For each 
raw separated image 7"B L 65MHz'( r )> we a PPly the convolution mask to the map to obtain a 
smooth plane for the scale of I = 2 1 pixels. Then, we apply the convolution mask to the 
obtained smooth plane to obtain a smooth plane for the next scale (i.e., / = 2 2 pixels). 
After repeating this process, a set of six smooth planes Ci(x,y) (i = 1, 2, . . . , 6) are 
obtained for scales of I = 2 n (n = 1, 2, . . . , 6) pixels, respectively. By calculating the 
difference between any two smooth planes that have adjacent scales, we obtain five subimages 
Wi(x,y) = Ci(x,y) — Ci + i(x,y), i = 1, 2, . . . , 5. The detection of clusters is performed on 
the subimages with the 3th and 4th scales, which are larger than the sizes of discrete sources 
convolved by the 21 CM A PSF. 

When the strength of the ICM magnetic field is set to be B = 2 fiG, 47 cluster-like 
structures are detected on the raw separated image 7"B L 65MHz'( r )5 °f which 37 are confirmed 
to be simulated clusters and 10 be simulated discrete sources, corresponding to a trust ratio of 
^detect /-^detect = 37/47 ~ 79%. We also have examined the central brightness temperature 
T CL 65MHz of the 37 detected clusters, and find that all of them have T CL 65MHz > 10 K, 
corresponding to a completeness ratio of A^ e L tect /Adulate I f^ L esMBa >w k = 37 / 48 - 77 % for a11 
simulated clusters that have Tq L 65MHz > 10 K (Fig. 7c). When B = 0.2 /iG is assumed for the 
ICM, 19 cluster-like structures (15 simulated clusters and 4 simulated discrete sources) are 
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detected, so that we have N$^/N™ = 15/19 ~ 79% and N^/N™^^^^ K = 
15/18 ~ 83%, respectively (Fig. 7(f). Discrete sources that are misidentified as galaxy 
clusters at this step can be recognized by examining their brightness temperature profiles 
and spectra (see details in §4.1). 

4. DISCUSSION 

4.1. Brightness Temperature and Spectral Information of Detected Bright 

Clusters 

In order to examine how accurate the FastICA + WT separation algorithm is, we 
randomly select three identified clusters with different redshifts and brightness temperatures 
(Table 2) to extract and study their radial brightness temperature profiles and spectra, and 
compare them with the input model, under the assumptions of B — 2 fiG and 0.2 /iG 
for the ICM, respectively. To correctly obtain a cluster's image in the 65 MHz channel 
(^B L 65MHz'( r ))5 nere we define a local circular region, which is large enough to cover the whole 
cluster (Fig. 7c & 7d), and feed the fake 21CMA survey maps extracted in this local region 
in five adjacent channels [y =61, 63, 65, 67 and 69 MHz, Av = 1 MHz for each channel) into 
the FastICA program. The separated images (Fig. 8) are obtained with the same FastICA 
+ WT algorithm and are rescaled in the same way as described in §3.2.1 . By repeating this 
process, we obtain separated images of each cluster in the 80, 95, 110, 125, 140, 155, 170, 
and 185 MHz [Av = 1 MHz for each channel) channels, respectively. Within a radius at 
which the brightness temperature drops to 1/10 of its central value, we extract the brightness 
temperature profiles and spectra of each cluster from the separated images (Fig. 9). At a 
given radius R as measured from the cluster's center ro, the systematic error is calculated 
as the mean deviation AT^J(R) = ^/E^I^V^ <f^J{R) >] 2 / N P ^\\ r - ro \ =R) 
where < T^ L J(R) > is the brightness temperature averaged at R, and -/Vpi xe i|| r _ ro |=R is the 
area. 

In Figure 9, we also show the original input radio brightness temperature profiles and 
spectra of the three selected clusters, as well as the difference between the extracted and input 
distributions as residuals. In each case we find that the extracted and input distributions 
agree very well with each other within the error bars. Consequently, we conclude that the 
FastICA + wavelet algorithm will enable us to separate and study the brightest (^og 5MHz > 
10 K) clusters in the 21CMA FOV, whose number is expected to be about 37 if B = 2 /iG 
(or 15 if B = 0.2 fiG). 

At 65 MHz, we successfully identify about 80% (37 out of 48 if B = 2 /iG, or 15 out of 
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18 if B = 0.2 /iG) of the simulated clusters that possess a central brightness temperature of 
T CL 65MHz = 10 — 5000 K, comparing to the average brightness temperature of ~ 2500 K of 
the foreground Galactic emission. The detected clusters are distributed in a mass range of 
2xlO 14 M -1.5xlO 15 M , andaredshiftrangeofz = 0.02-1.6 if B = 2 fj,G, or z = 0.02-0.8 
if B = 0.2 fiG. All of them are bright enough for a direct model study to determine their 
spectral properties. We ascribe the failure of detection of the remaining 14 clusters (11 
clusters if B = 2 /iG, or three clusters if B = 0.2 /iG) to the following three reasons. First, 
six out of the 14 undetected clusters exhibit relatively flat spectra, which makes it difficult 
to disentangle them from the Galaxy component, the spectral index of which ranges from 
~ 2.7 to 2.9. Second, five out of the 14 clusters are faint and less extended, either embeded 
in brighter clusters or overlapped by bright discrete sources. And third, three out of the 14 
clusters are less extended and are masked along with adjacent bright discrete sources. We 
also find that 10 (if B = 2 /iG) or four discrete sources (if B = 0.2 fiG) are misidentified 
as galaxy clusters. This mistake can be corrected by examining the brightness temperature 
profiles and spectra of the sources. 

For the cosmological 21cm signals survey projects such as the 21CMA, although our 
separation approach can only extract bright clusters, which are more than two orders of 
magnitude brighter than the cosmological 21cm signals, from the images, the analysis on 
the cluster component can deepen our understanding on cluster characters (e.g., brightness, 
shape, spatial distribution, spectrum) in the low-frequency radio sky, the information of 
which is valuable when detecting the cosmological 21cm signals in both image and frequency 
spaces. 



4.2. Existence Rate of Cluster Radio Halos 

In the frame of the most popular model for cluster radio halos, i.e., the turbulent re- 
acceleration model (e.g., Ferrari et al. 2008), every cluster should host a radio halo after 
experiencing a major merger. However, since the halo's radio emission decays quickly as the 
relativistic electrons lose their energy, radio halos are expected in only about 20 — 30% of 
massive (M CL > 2 x 10 15 M ) clusters and in 2 - 5% of less massive (M CL ~ 10 15 M ) 
clusters (Cassano et al. 2006), in other words, in the whole sky < 50 clusters are found to 
host a bright radio halo at 1.4 GHz according to the obsevational works (Cassano et al. 2006; 
Giovannini et al. 2009; Brunetti et al. 2009). Using this halo number density we estimate 
that < 1 radio halo shall appear at 1.4 GHz in the 21CMA FOV, which is consistent with the 
1.4 GHz detection of Abell 2294 in this field. On the other hand, by calculating the evolution 
of the power density ghz(^e) °f cluster radio halos (see Appendix), we find that among 
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1084 simulated clusters only one (if B = 2 fiG) or none (if B = 0.2 fiG) cluster (Fig. 10) 
is brighter than the standard VLA detection threshold at 1.4 GHz, given an intergration 
time of ~ 3 hrs (Giovannini et al. 2009). This shows that our work is nicely consistent with 
current observations performed at higher frequencies, and the prediction for the apperence 
of more than one dozen low-frequency cluster radio halos in the 21CMA FOV is reliable. 

This work is accomplished based on the predictions of the standard turbulent re- 
acceleration model, which assumes that most cluster radio halos are formed in major merger 
events. However, in alternative models (e.g., the magneto- hydrodynamic (MHD) model; 
Brunetti et al. 2004; Cassano 2009), some radio halos, although relatively small and weak, 
can also be produced on cluster scales via minor mergers, so that much more halos are ex- 
pected in the low-frequency radio sky. The survey of galaxy clusters with upcoming facilities 
working in the band of a few hundred of MHz will help distinguish these models. 

5. SUMMARY 

By carrying out Monte-Carlo simulations we model the foreground of the cosmological 
reionization signals, including emissions from our Galaxy, galaxy clusters, and extragalactic 
discrete sources (i.e., star- forming galaxies, radio-quiet AGNs, and radio-loud AGNs). As an 
improvement of previous works, we consider in detail not only random variations of morpho- 
logical and spectroscopic parameters within the ranges allowed by multi-band observations, 
but also evolution of radio halos in galaxy clusters, assuming that relativistic electrons are 
re-accelerated in the ICM in merger events and lose energy via both synchrotron emission 
and inverse Compton scattering with CMB photons. With these we create the fake 50 — 200 
MHz sky as will be observed with the 21CMA array. We show that about 80% (37 out of 
48 if B = 2 fiG, or 15 out of 18 if B = 0.2 /zG) of the clusters with a central brightness 
temperature of > 10 K at 65 MHz can be successfully disentangled from the fake sky maps, 
if an ICA+wavelet detection algorithm is applied. 

We thank the referee for the valuable comments on our manuscript. This work was 
supported by the National Science Foundation of China (Grant No. 10673008, 10878001 
and 10973010), the Ministry of Science and Technology of China (Grant No. 2009CB824900, 
2009CB824904), and the Ministry of Education of China (the NCET Program). 
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Appendix: Evolution of Synchrotron Emission Spectra of Cluster Radio Halos 

By using the approach presented in You (1998; p201), we calculate the evolution of 
synchrotron emission spectrum for a population of relativistic electrons that lose energy via 
synchrotron radiation and inverse Compton scattering with CMB photons. Assuming that 
a group of relativistic electrons are injected into the magnetic field at the initial time (i.e., 
t — 0) once only, the number of which possesses an energy distribution of 

N( 1 ,t)\ t =o = N % r , (20) 

where N is the density normalization, 7 is the Lorentz factor so that the electron energy 
is 7m e c 2 , 70 = j\t=o is the initial electron energy, and T (defined to be 3 in this work) 
is the initial energy spectral index. Due to the conservation of electron number, electrons 
distributed within the energy interval 7o ~ 7o + dlo at £ = shall fall in the energy interval 
7 ~ 7 + g?7 at a later time time t, we have 

N( 1 ,t)d 1 = N % r d l0 , (21) 

where iV(7, t) is the electron number per unit volume at time t within the energy interval 
7 ~ 7 + drf. We also know that electrons lose energy via synchrotron radiation and inverse 
Compton scattering with CMB photons at a rate of 

7 = -6 7 2 , (22) 

where b = 3 x 10 8 £/ mag , and C/ mag = [B 2 + (B CMB ) 2 ]/8ir (B is the local magnetic field 
strength of the ICM, and _B CMB is the equivalent magnetic field strength due to the inverse 
Compton scattering between relativistic electrons and CMB photons) is the energy density 
of the synthesized magnetic field (see also Parma et al. 1999; Ferrari et al. 2008), by solving 
which we obtain 

7 = 7 ° , (23) 
1 1 + &70* 

or 

To = (24) 

and 

d7 ° = (i^W d7 " (25) 

Then by substituting equation (24) and equation (25) into equation (21), we have 

iV(7, t) dl = I iVo7 ~ r(1 - 67t)r " 2rf7 ' 71 < 7 < 72 ' (26) 
^ ' I 0, 7 < 7i or 7 > 72 , K ' 
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where 71 and 72 are the lower and upper limits of electron energy at time t, respectively. 

On the other hand, if the synchrotron radiation of the electron population is isotropic, 
the synchrotron emission coefficient (the energy emitted per unit time per unit volume per 
unit frequency) can be simply expressed as 

>(*)=4tt J N(n,t)PZdri, (27) 

where Pj; oc f^J„ K 5 / 3 (x)dx is the power density of one electron's emission, = is 
the Larmor frequency, v c is the critical frequency, and K 5 / 3 is the modified Bessel function of 
order 5/3. Therefore, the power density (the energy emitted per unit time per unit frequency) 
at any frequency v and time t can be calculated with both equation (26) and equation (27), 

P?\t) = Ut)V = m ^^l = °l (28) 

JlA GHz(r = 0) 

where V is the volume of the cluster radio halo, and the initial power density at 1.4 GHz, 
i.e., P X C 4 GHz (t = 0), has been derived from the X-ray temperature of the cluster (§2.2.2). 
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Table 1: Estimated numbers of extragalactic discrete sources and their structural subcom- 
ponents in the 21CMA FOV (§2.3) 



Source type 


Number of sources (10 6 ) 


Number of Structures (10 6 ) 


Radio-quiet AGN 


7.1 


7.1 


FR I 


4.7 


14.1 


FR II 


0.0005 


0.0025 


Normal galaxies 


40.8 


40.8 


Starburst galaxies 


1.4 


1.4 


GPS & CSS 


0.4 


0.4 



Table 2: Confusion level estimation with different methods (§2.4)t 



Method 






85 MHz 




1. 


50 MHz 










21CMA Instrument 


& Our Image (PSF= 


=3.2') 


21CMA Instrument (PSF = 


=1.6') 


Our Im 


age (PSF=3.2') 


I 


Our flux density function 
Literature flux density function 


36.1 
33.1 


m.Iy/407 K 
mJy/373 K 




3.33 mJy/28.2 K 
3.20 mjy/27.1 K 




18.4 
18.2 


mJy/39.2 K 
mJy/38.8 K 


II 


Photometric Criterion 


32.0 


mJy/361 K 




3.65 mJy/30.9 K 




18.6 


mJy/39.6 K 


III 


SExtractor 


28.1 


mJy/317 K 




3.02 m.Iy/25.6 K 




14.1 


mJy/30.0 K 



The confusion level given are 3<r values, corresponding values of brightness temperature are also listed. 



Table 3: Properties of the randomly selected three clusters, which are separated from the 
fake 21CMA field (§4.1)* 



Cluster ID 


M CL 


z 


^200 




A C 4GHz(* = 0) 


t E 


rfCh 

1 0, 65MHz 




(10 14 M ) 




(ftfj 1 Mpc) 


(keV) 


(10 23 W Hz" 1 ) 


(Gyr) 


(K) 
















B = 2 fiG B = 0.2 fiG 


1 


3.02 


0.07 


1.35 


3.7 


4.37 


0.24 


636.0 129.2 


2 


3.04 


0.25 


1.27 


3.9 


5.09 


0.11 


146.0 32.5 


3 


8.49 


0.27 


1.78 


7.0 


29.0 


0.11 


634.1 111.9 



* Cluster identifications, masses, redshifts, virial radii, X-ray temperatures, initial power densities at 1.4 
GHz, evolving time, and central brightness temperatures at 65 MHz (if B = 2 ^G or 0.2 fiG) arc listed, 
respectively. 
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Fig. 1. — Simulated two-dimensional images for the distributions of (a) brightness tempera- 
ture of the Galactic synchrotron emission at 65 MHz, (b) Galactic synchrotron temperature 
spectral index (a^ syn ) , and (c) brightness temperature of the Galactic free-free emission at 
65 MHz, all constrained in the 21CMA FOV. 




Fig. 2. — The theoretical number density distribution n(M CL > M, z) (count per unit 
redshift per unit solid angle) of galaxy clusters as a function of mass and redshift according to 
the Press-Schechter model. In our simulation, we define n CL (z) = n(M CL > 2 x 10 14 M , z). 
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Frequency (MHz) X-ray temperature (keV) 



Fig. 3. — (a) Spectral evolution of cluster radio halos at z = 0.2, assuming B = 2 /iG and 
B = 0.2 fiG, respectively, (b) An approximation to the relation between PfJ ghz(^ = 0) an d 
Tx (solid line), which is obtained by repeatedly fitting the observed data with maximum 
likelihood method and excluding the data located outside the 3a lower limit ($2.2.2). The 
marks "o"s and "x"s stand for the data finally involved and not involved in the fitting, 
respectively. We also superpose the data of three new cluster samples provided by Giovannini 
et al. (2009) as a comparison (black triangles). 
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Fig. 4. — Cluster brightness temperature images simulated at 65 MHz, assuming B = 2 /zG 
(a) and B = 0.2 //G (b), respectively. Both of the images have been convolved with a PSF 
of 3.2' and are constrained in the 21CMA FOV. 
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Fig. 5. — Brightness temperature images of extragalactic discrete sources simulated at 65 
MHz, which have been convolved with a PSF of 3.2', and are constrained in 2° x 2° fields to 
illustrate structural details. 
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Fig. 6. — The total flux density functions of our simulated extragalactic discrete sources. 




Fig. 7.— (a)-(b): Simulated 65 MHz 21CMA maps, assuming B = 2 jiG and B = 0.2 
fiG, respectively, (c)-(d): Cluster component separated at 65 MHz, where the green squares 
mark the correctly detected bright (Tqq 5 MHz > 10 K) clusters, the orange squares mark 
the missed bright clusters, and the red squares mark the sources misidentified as clusters 
(§3.2.2). Three magenta circles mark the clusters randomly selected to examine their spatial 
and spectral properties (§4.1). 
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Fig. 8. — Comparison between the simulated (left panels) and separated (right panels) 
65 MHz brightness temperature images of the three randomly selected clusters, assuming 
B = 2 fiG and B — 0.2 fiG, respectively. The green circles mark the radii at which the 
brightness temperature drops to 1/10 of the central value. 
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Fig. 9.— Brightness temperature profiles in 65 MHz channel and spectra of the three 
randomly selected clusters extracted within the green circles defined in Fig. 8 (crosses), as 
compared with the simulated results (solid lines). Residuals obtained in the \ 2 test between 
the separated distribution and simulated distribution in each case are also plotted in the 
lower panels. 
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Fig. 10. — -Pi°4 ghz(^) — z distribution of the simulated clusters, assuming B = 2 fiG ("*"s) 
and B = 0.2 fi G ("<}"s), respectively, as compared with the VLA sensitivity limit, assuming 
a standard observation with an integration time of ~ 3 hrs (Giovannini et al. 2009). The 
star marks Abell 2294, which is the only one detected cluster that holds a radio halo in the 
21 CM A FOV nowadays. 



